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Abstract 

In this note we propose an overview on the current theoretical and experimental limits on a Higgs singlet 
extension of the Standard Model. We assume that the Boson which has recently been measured by the 
LHC experiments is the lightest Higgs boson of such model, while for the second Higgs Boson we consider 
a mass range of 600 GcV < win < 1 TeV, i.e. outside the range of the direct searches presented so far. In 
this light, we study the impact of perturbative unitarity limits, Renormalisation Group Equations analysis 
and experimental constraints (Electro- Weak Precision Tests, measurements of the light Higgs coupling at 
the Large Hadron Collider). We show that, in the case of no additional hidden sector contributions, the 
largest constraints for higher Higgs masses stem from the assumption of perturbativity as well as vacuum 
stability for scales on the order of the SM metastability scale, and that the allowed mixing range is severely 
restricted. We discuss implications for current LHC searches in the singlet extension, especially the expected 
suppression factors for SM-like decays of the heavy Higgs. We present these results in terms of a global 
scaling factor n as well as the total width L of the new scalar. 
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I. INTRODUCTION 

The recent discovery of a new particle [H [2] which is in accordance with the scalar Boson from the 
Higgs mechanism [3HZ] by the LHC experiments is one of the big breakthroughs in contemporary 
particle physics. If the discovered particle is indeed the Higgs Boson predicted from a SM-like 
Higgs-doublet sector, all its properties are completely determined by theory. Therefore, the 
current quest of the theoretical and experimental community is to establish whether the properties 
of such particle are in accordance with standard predictions, or it is only a component of a more 
involved Higgs sector. For this, all couplings as well as the spin structure of the new particle need 
to be severely tested. 

In this work, we consider the simplest extension of the SM Higgs sector, i.e. we add an 
additional singlet which is neutral under all quantum number of the SM gauge groups [U [9] and 
acquires a vacuum expectation value ( VEV) . This framework has already been discussed by several 
authors |10H22j . and can easily be interpreted in terms of one of the benchmark models of the 
heavy Higgs cross-section working group [23J. In this work, we investigate its extended parameter 
space with the assumption that the heavy Higgs mass lies in the range 600 GeV < ran < 1 TeV 
and there are no relevant contributions from hidden sectors. We combine the effects of several 
constraints: LHC bounds on the light Higgs signal strength, bounds from perturbative unitarity, 
electroweak (EW) parameters in terms of S,T, and U, and limits from perturbative running of 
the couplings. As a major result, we find that for mn ^ 700 GeV especially the running of the 
couplings severely restricts the allowed parameter space of the model, leading to scaling factors 
in the percent range. We comment on the effects of these restrictions on LHC production cross 
sections at a 8 and 14 TeV LHC as well as on the SM-like branching ratios of the model. 

This paper is organized as follows: In Section II, we briefly review the model setup. Section III 
is devoted to the investigation of the allowed parameter space taking all constraints into account. 
In Section IV, we comment on the impact of these limits on LHC observables. We summarize in 
Section V. 
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II. THE MODEL 



A. Potential and couplings 

In this paragraph we will shortly review our model: we enlarge the SM Higgs sector with a 
further real Higgs singlet x-, which is pure singlet under each gauge group of the SM [8j 19} 124"]. 
The most general gauge-invariant and renormalisable scalar Lagrangian is then: 

% a = (D"fO t D » H + (^'X) t D^ X ~ V(H, x) , (1) 
with the scalar potential given by 

V(H, x) = -m 2 H^H -/i 2 | x I 2 + ( H^H 




-m 2 



H^H - n 2 | x I 2 +Ai(# t #) 2 + A 2 | X I 4 +\^H | x | 2 , (2) 



where x is the Vacuum Expectation Value VEV associated to the new Higgs field. We here 
implicitely impose a Z2 symmetry which forbids additional terms in the potential. 

To determine the condition for V(H, x) to be bounded from below, it is sufficient to study its 
behaviour for large field values, controlled by the matrix in the first line of Eqn ([2]). Requiring 
such a matrix to be positive-definite gives the conditions 

4AiA 2 -A^ > 0, (3) 
Ai,A 2 > 0, (4) 

where the condition given by Eqn. Q corresponds to the requirement that the potential is bounded 
from below for large field values, while Eqn. ([3]) guarantees that the extremum is indeed a local 
minimum. 1 Since the physical mass eigenvalues are gauge invariant, we define the Higgs fields 
following the unitary-gauge prescription: 

I \ 



H = 



h + v 



_h' + x 



1 We give the exact derivation of the resulting eigenstates and the derivation in Appendix [X] and here only cite the 
relevant results. 
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The explicit expressions for the scalar mass eigenvalues are: 

mi = \ lV i + A 2 x 2 - v 7 ^ 2 - A 2 x 2 ) 2 + (X 3 xvy, (6) 
m 2 H = Xiv 2 + X 2 x 2 + ^{Xiv 2 - A 2 x 2 ) 2 + {X 3 xv) 2 , (7) 

where h and H are the scalar fields of definite masses and ran respectively, with m\ < m 2 H . 
These eigenvalues are related to the following eigenvectors: 

cos a — sin a \ I h 

sin a cos a 

where -f < a < § fulfils 2 : 




(8) 



sin 2a = — , 9 

V(Ai^ 2 - A 2 x 2 ) 2 + (X 3 xv) 2 

o - X x v 2 

cosza = — . (10) 
V(Ai« 2 - A 2 x 2 ) 2 + (A3:™) 2 

From Eqn. ([8]), it is clear that the light (heavy) Higgs couplings to SM particles are now suppressed 
by cos a (sin a). 

From equations @-0-([9]), it is straightforward to have: 

A i = ttI + m fe) sin 2 a _ ^ cos 2 a + sin 2 a 

2u z 2v z 2v z 2ir 

A2 = 2^ + 2x 2 COS a= 2^ Sm a+ 2# COS a 

a 3 = 4^-.w (id 

In summary, the heavy Higgs is a "twin" version of the light Higgs with rescaled couplings to 
the matter contents of the SM. In fact, the only novel channel with respect to the light Higgs case 
is H — > hh. The decay width T and coupling strength a! of the H — >• hh decay are [U I24j : 

8 it ran V m z H 

y! = — ^ (x cos 3 a + v sin 3 a) + (A3 — 3 Ai) v cos 2 a sin a + (A3 — 3 A 2 ) x cos a sin 2 a. 

(12) 



We here briefly discuss the behaviour of when x and sin a are varied: from Eqn. (12), it 
is clear that 3 |/x'| 2 (| sin a|) > |//| 2 (— sin a). The difference is more pronounced as x is increased. 



In all generality, the whole interval < a < 2n is halved because an orthogonal transformation is invariant under 
a — > a + ty. 

3 If x 3> v, we can approximate /j,' m A3 2; cos a (sin 2 a — \ cos 2 a) — ?>\\v cos 2 a sin a + O (^). We then have 
//(sin a) = — fi'(— sina). If the terms ~ ©(a; -1 ) cannot be neglected, they introduce a positive/ negative 
contribution to \n'\ depending on the sign of sina. 
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In addition, for a fixed value of sin a > 0(sina < 0), |//| decreases (increases) constantly for 
increasing x. These features will become important in the discussion of the experimental and 
theoretical constraints in the next sections. 

B. Number of free parameters 

Our simple singlet extension model has in principle 5 free parameters on the Lagrangian level 

Ai, A2, A3, v, x. 

The coupling parameters Aj are related to the masses and the effective mixing according to Eqns. 
(|6]),([7]), and Q, and we obtain the independent parameters 

m-h, mH> o, v, x. (13) 

Moreover, we will reexpress x by tan f3 according to 

n V 

tanp = — 

x 

to accomodate for standard notation in models with exended Higgs sectors. If we assume the 
vacuum expectation value of the Higgs doublet value to be Standard Model-like such that v ~ 
246 GeV, and equally set the Higgs mass of the light Higgs to = 125 GeV, we are left with 
three independent parameters itlh, a, (x/ tan/3). All results in the following sections will be given 
in dependence on these variables. In this work, we restrict the range of the Singlet VEV to 
x G [100 GeV; 10 TeV] , leading to tan G [0.025; 2.46]. 

III. THEORETICAL AND EXPERIMENTAL BOUNDS ON THE HIGGS SINGLET EX- 
TENSION 

In this section, we will discuss the current theoretical and experimental limits on the singlet ex- 
tension model. We here consider 

• limits from perturbative unitarity, 

• limits from EW precision data in form of the S, T, U parameters, 

• perturbativity constraints on the couplings, as well as conditions on a potential which is 
bounded from below, 
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• limits from measurements of the light Higgs signal strength, 

• limits from perturbativity of the couplings as well as vacuum stability up to a certain scale 
/i mn , where we chose /x run ~ 10 10 GeV, 10 19 GeV as benchmark points. 

In this chapter, we will investigate the parameter space (sin a, tan/3), while keeping mn fixed; 
however, in order to demonstrate the effects of the partial-wave treatment of perturbative unitarity, 
we will equally comment on the highest possible mass of the heavy Higgs mjj max m this parameter 
space, including exclusion bounds from electroweak precision data using m^max- We discuss all 
limits separately in the following subsections. 



A. Limits from perturbative unitarity 

Tree-level perturbative unitarity |25] [26] puts a constraint the Higgs masses of our theory via a 
relation on the partial wave amplitudes aj(s) of all possible 2 — > 2 scattering processes: 



|Re( aj ( S ))| < \. (14) 

In the high energy limit, \fs — > oo, only the ao partial wave amplitude does not vanish, 
instead it approaches a value depending only on m^, ttlh, a and x. Therefore, by applying the 



condition in eq. (14), we can obtain several different (correlated) constraints on the Higgs masses 
and mixing angle, i.e., we can find the m^-m^-a subspace in which the perturbative unitarity 
of the theory is valid up to any energy scale. We therefore studied the unitarity constraints in 
our model by calculating tree- level amplitudes for all two-to-two processes 4 X\ X2 — > Y1Y2, with 
(X\,X2), (11,12) 6 (W + W~ , ZZ,hh,hH,HH) in terms of the mixing angle between the two 
physical Higgs fields and their masses. Then, we calculated the normalized 5-dimensional bosonic 



scattering matrix and we imposed the condition of Eqn. ( 14 ) to each of its eigenvalues (the largest 
in modulus gives the best constraint). Note that, in accordance with |24j . the constraint based 
on generic unitarity considerations (cf eg |15] ) for the heavy Higgs of m# < 700 GeV is much 
loosened 5 . 



4 Calculations where actually carried out with the vectors Bosons being replaced by the corresponding Goldstone 
Bosons following the equivalence theorem [27] . 

5 This result is also confirmed in [28], where a similar scenario is investigated. 
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Maximally allowed heavy Hlggs masses from perturbative unitarity Maximally allowed heavy Higgs masses from perturbative unitarity 




sin a sin a 



FIG. 1: Allowed mixing ranges for different maximal heavy Higgs masses m/f ;max in dependence of the 
mixing angle a as well as ratio of the VEVs tan/3. LEFT: Heavy Higgs masses above (1, 2, 3, 4) TeV are 
excluded in the (red, green, blue, magenta, yellow) area (from sin a = —1 to sin a = 0). The region 
exluded for small sin a and tan/3 > 2 leads to mjj imax < 600 GeV. The most important constraints on 
the upper limit in the small mixing region stem from scattering processes involving heavy Higgs pairs (not 
shown here). RIGHT: Zoom into the region where TOH,max < ITeV. 

Figure [T] shows the regions in parameter space which are still allowed after limits from per- 
turbative unitarity only. We found that for small mixing angles within our scan range, the most 
dominant contribution stems from scattering processes involving only heavy Higgses. For sin a ~ 0, 
the scattering matrix becomes approximately block diagonal with a SM block and the decoupled 
HH — > HH element, and the latter gives the unitarity limits on the singlet VEV, i.e. 6 

16 7T I! 2 

tan 2 /3 < =- + 0(a) for a (HH -> HH) < 0.5; 

if tan f3 is decreased accordingly, this boundary can therefore be fulfilled any heavy Higgs mass. 
For small, but non-zero mixing angles |sina| ~ 0.02 and tan/3 < 0.1 , upper limits for the 
allowed maximal heavy Higgs mass can reach up to 35 TeV. We found that generically, heavy Higgs 
scattering processes dominate for tan (3 > 1.5; if tan j3 is decreased, and for non-zero mixing, gauge 
Boson scattering becomes equally important. However, in most cases the whole 5x5 scattering 
matrix involving all partial wave contributions needs to be considered, and an approximation 
considering a single dominant process cannot give a valid predicition of the upper limit on the 
allowed heavy Higgs mass. 

6 This boundary is in fact stronger than perturbativity of the coupling alone, which leads to tan 2 ft < 8 ^ for 
sin a = 0. 
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B. Limits from electroweak precision data 

Constraints from EW precision data are incorporated using the S, T, U parameters (29] EQ] , 
which parametrize deviations from the SM predictions and thereby render constraints on new 
physics from higher order corrections stemming from BSM contributions. 

We here follow [31], which cites the values for the EW parameters as 

S = 0.00 ± 0.10, T = 0.02 ± 0.11, U = 0.03 ± 0.09 

and equally used = 125 GeV as an input value for the calculation of the SM reference values 
where S ie f = T re f = U rc { = 0. As a cross check, we have compared results from our code with 
the values for S, T, U for all benchmark points specified in |24j . and found agreement with small 
variations on the 10% level (we used 7 (rn% op ,mh) = (173.5, 150) GeV). To accomodate for this 
slight disagreement, we decreased the allowed regions for S, T, U to 

S = 0.00 ± 0.095, T = 0.02 ± 0.105, U = 0.03 ± 0.085 

in our scans, and use as input variables |32j 

s z = 0.2313, a s {M z ) = 0.120, m t = 173 GeV. 

We then use [33] 

-Xtot = cos 2 ^^^^) + sin 2 a X(i71h) 

with X € [S, T, U] . Note that this approach neglects suppression of the couplings in all 
but the leading order, and equally does not take the Hhh couplings into account which can 
appear in higher order corrections including heavy Higgses running in the loops 8 . However, 
as we will argue below, EW precision data basically poses no constraint on the parameter 
space after all other restrictions have been taken into account. Of course, a more detailed 
analysis would be desireable here, and is in the line of future work. In our approximation, 
the constraints basically rule out values of |sina| > 0.5 — 0.7 depending on tan/3, where the 
strongest constraints here come from the T-parameter. U does not pose any additional constraints. 



7 In [21], the actual value of the top mass which was used is not given. Variations for m t £ [170.5; 173.5] did not 
significantly change our results. 

8 Cf also [12] for a generic calculation with multiple scalar extensions of the SM. 
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This closes our discussion of scans using maximally allowed heavy Higgs masses from 
perturbative unitarity. We found that within our scan range the maximally allowed 
Higgs masses are 0(35 TeV) for small mixing angles and that, using these maximal Higgs 
masses, EW precision data give additional constraints in the large mixing regions. For 
mH,max < lTeV(2TeV), |sina| < 0.6(0.5). For m# ima x > 2TeV, EW precision data give no 
additional constraints to our model. 

In the following, we will fix the Higgs mass to mn £ [600, GeV; 1 TeV], and equally include the 
measurement of the light Higgs signal strength \u\ as well as vacuum stability and perturbativity 
of the couplings up to a metastable scale/ the Planck scale. We will see that indeed these latter 
requirements are much more stringent than EW precision data and the light Higgs measurements 
and render severe constraints on the parameter space of our model. 

C. Constraints from the signal strength of the light Higgs 

If we want to accomodate for the light Higgs measurements [TJ |2|, we need to take into account 
the limits on the maximally allowed value of | sina| from the overall signal strength \u\. In general, 
we have 9 

= cr B SM(m h ) _ cos 4 aT S M(m h ) 

csM^ft) cos 2 ar SM (m),) + sin 2 aT hid (m h ) ' 

where r SM / hid( m h) denote the decay widths of the light Higgs in the SM/ a possible hidden 
sector. We here constrain ourselves 10 to cases where T^id = 0; then, the above equation leads to 
[i = cos 2 a. 

The values measured by the LHC experiments [TJ [2] then render 11 

|sina| e [0;0.23] from /i G [0.95; 1]. 

9 In fact, loop-induced couplings like the h — ¥ 77 branching ratio in principle call for a more refined treatment, 
cf eg discussion in |21j . However, the corresponding corrections are generally on the sub-permill level, and can 
therefore safely be ignored in our simple limit-setting. For fitting procedures, on the other hand, such a more 
complex coupling structure needs to be taken into account. 

10 Taking decays in the hidden sector into account additionally reduces the allowed mixing range. To understand 
this, consider the case that Tsm = rW; in this case, the constraint is strengthened to cos 4 a > 0.95, leading to 
sin a < 0.17. The case in which Thid = is therefore the best case scenario. 

11 The official ATLAS fit for a 126 GeV Higgs are given by [TJ [i = 1.4 ± 0.3, and from CMS for a 125 GeV Higgs 
as [2] fj, = 0.87 ± 0.23. Being conservative, we consider fi > 0.95, where we take the fact into account that the 
model considered here cannot accomodate for fj, > 1. We also assumed that the errors of the ATLAS and CMS 
measurements are completely uncorrelated. Newer values 34, 35 do not significantly change this result. Cf also 
|21| for a best fit result for this model. 
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Effects of perturbative unitarity limits for m H = 600 GeV Effects on parameter space from PU varying the Higgs mass 




sin a sin a 



FIG. 2: LEFT: Allowed regions in parameter space where mn = 600 GeV. The only restriction comes 
from PU, which gives an upper limit on tan/?. RIGHT: Allowed parameter space for a 600 GeV as well as 
1 TeV heavy Higgs. As before, the only limit in parameter space comes from perturbative unitarity, which 
fixes the upper allowed bound of tan j3 to 1.9 and 1.0 respectively. 

These limits on the measurement of the 125/126 GeV Higgs Bosons coupling strength are in fact 
much more stringent than EW precision observables. 

Figure [2] shows the limits for ran = 600 GeV and ran = 1 TeV respectively. The only constraint 
arises here from perturbative unitarity, which sets an upper limit on tan/3 in both cases. This is 
generically due to a large A2 value in these regions of parameter space. In accordance with the 
behaviour of A2 ~ ra\ for fixed (sin a, tan/3) values, we equally observe that the coupling gets 
larger for larger ran values, leading to a decrease in the upper limit of tan /3. Most of the parameter 
space ruled out by perturbative unitarity would however be equally excluded by the requirement 
of perturbativity of A2 at the EW scale, as discussed below. 



D. Limits from perturbativity and vacuum stability 

We equally consider vacuum stability as well as perturbativity of the Higgs potential couplings 
up to a certain scale /z mn . Vacuum stability follows from Eqn. Q, while perturbativity of the 
couplings leads to the requirement that 

Al,2(/^mn) < 4-7T, |A 3 (/i run )| < 47T. 

At the electroweak scale, we found that these conditions pose no additional constraints on the 
allowed parameter space of the model, when limits from the light Higgs signal strength and per- 
turbative unitarity are taken into account. If we neglect perturbative unitarity limits, the upper 
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allowed values of tan/3 following from perturbativity of the couplings alone are 2.05(1.24) for 
mji = 600GeV(lTeV) (for sin a = 0) respectively, which slightly enhances the allowed tan/3 
ranges. Before considering the running of the couplings, we can therefore say that 

• perturbative unitarity alone indeed allows for heavy Higgses in the 30 TeV range 

• the strongest constraints considered so far, when the experimental results for the light Higgs 
signal strength are taken into account, stem from perturbative unitarity. 

We now discuss limits from perturbativity up to /i run , where we use the running parameter 
t = In ^^f 1 ^ such that t = for /x run = v. We here impose the constraint given by Eqn. (3) at 
all energies. Note that in a strict sense this is not required for vacuum stability; for positive A3 
values, fulfilling Eqn. Q is sufficient, c.f. e.g. the discussions in [36j|37]. However, as we require 
perturbative unitarity up to arbitrary high scales, we also demand that the process of electroweak 
symmetry breaking remains the same and that therefore the minimum of the potential is indeed 
positioned at the VEVs of the two fields; this approach has e.g. been followed in [38]. We will 
briefly comment on the effects of releasing such a condition on the collider observables in Section IV. 



In the following discussion, we mostly focus on mu = 600 GeV, t = 37, but will equally 
give results for mu = 1 TeV and /i run = 10 19 GeV. In the end of the discussion, we will 
comment on the generic changes for a higher Higgs mass or the requirement of perturbativity 
and vacuum stability at higher scales. For the sake of the argument, we will temporarily neglect 
the measurement of the light Higgs signal strength and consider mixing angles | sina| < 0.49 in 
the discussion of RGE running effects at the low scale, in order to exemplify the generic effects 
on the parameter space. The signal strength measurement will however be included again in the 
discussion of collider observables in Section IV. 



The renormalization group equations for this model are given by |24l [38] 
1 



d 
dt 
d 

dt 
d 

dt 



Ai 

A 2 
A3 



16vr 2 1 2 

1 



1 A 2 + 12 A? + 6 A! yl - 3 yf - \ X, (3 5 2 + g\) + A [2 g* + ( 5 2 + gf) 2 



16 7T 2 
1 

16 vr 2 



[A3 + 10 A 2 ,] , 



A, 



6A 1 + 4A 2 + 2A 3 + 3y 2 -^(3r 7 2 + 9 2 ) 



where y± is the (equally running) top Yukawa-coupling and g, g\ are the running couplings of the 
SM gauge groups. For the decoupling case as well as for cross check for the running of the gauge 
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couplings, for which we chose the analytic solution at one loop, we reproduced the results in [39] . 
where the SM breakdown scale following the one-loop treatment here was at t = 36 corresponding 
to a scale /i mn ~ 1.6 x 10 10 GeV. 

We found that the strongest constraints from a phenomenological viewpoint, ie upper limits 
on the allowed mixing angle, actually stem from perturbativity of Higgs self-couplings Ax, A2; for 
t = 37 and low mn, we found the requirement that | sina| < 0.3. This poses a much stronger 
constraint than electroweak precision tests. In the following, we discuss limits from perturbativity 
as well as vacuum stability in more detail: 

• Perturbativity of Ai and upper limit on | sina| 

The strongest constraint on large mixing angles stems from the running of Ai . For tan /3 < 
0.1 and large mixing angles 



at the electroweak scale, so the /3 function of Ai is positive. In this case, Ai quickly grows 
and approaches the upper limit of 4tt (eg for sin a ~ 0.49, tan/3 ~ 0.025, this is reached 
for the relatively low scale of /i run ~ 350 TeV). This remains the dominant effect until 
tan/3 > 0.36, where A2 starts to rise more quickly. From the running of Ai, we obtain 
I sina| < 0.3 (0.2) for running up to t = 37 (the Planck scale) for a 600 GeV Higgs; for 1 
TeV, these values change to 0.2 (0.12). 

• Perturbativity of A2 and upper limit on tan /3 

For tan /3 > 0.36, the most dominant constraint comes from the running of A2 in almost all 
regions of parameter space 12 . Generically, a good estimate of the limits can be obtained by 
considering the zero-mixing case and tan/3 > 1: we then have 



In this case, it is easy to estimate the maximal value of tan/3 allowed such that A2 = 4n. 
The corresponding /3-function can be reduced to 



A 2 , I A3 1 <^ Ai 



A2 ~> Ax, A3. 



d\2 
dt 




which has the solution 



A 2 (t = 0) 



1 - 



^t\ 2 (t = 0) 



12 



In the region where | sina| > 0.26, A3 running sets in as well, cf discussion below. 
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Requiring \2(t) < A max then leads to 

Amax — A(t = 0) 



t < 



f^Mt = 0) A maj 
As A 2 (i = 0) ~ this translates to a lower limit on x 



2 

2 { \ H 



1 5 

1 + 

Amax O 7T 



(here we set sin a = 0) and therefore tan/3 max = — — . Inserting explicit values for mu = 

v ^min 

(600, GeV; 1 TeV) for t = 37(76) renders (tan/3) max = (0.37(0.26); 0.22(0.15)). These 
values agree with our numerical findings. 

Perturbativity of A3 and restriction in the large tan/3/ large sin a region 
In a small region for tan j3 ~ 0.4 and large positive mixings, 

A3 ^ Ai, A 2 , 

which corresponds to the transition between Ai and A2 dominance. In this region, all cou- 
plings evolve similarly fast up to high scales. As an example, we show the running of all 
Higgs sector as well as the top Yukawa coupling for a point in this part of parameter space 
in Figure [3j 

First vacuum stability condition fAi > 0) and minimal mixing angle \ sinctj 
For small (or 0) mixings, this is the well-known metastability problem of the SM Higgs 13 with 
a low mass of 125 GeV. In our scan, the couplings becomes negative at a scale t = 36, which 
corresponds to roughly o Tnn ~ 1.6 x 10 10 GeV. For small mixing angles |sina| < 0.001, 
the problem persists. There is no significant change from this limit for raising the Higgs 
mass to 1 TeV 14 . 

Third vaccum stability condition (4Ai A2 > A|J and minimal mixing angle \ sina| 
For 0.001 < I sina| < 0.4 and tan/3 < 0.4, the third vacuum stability condition 

4 Ai A 2 - A3 > 

poses the largest constraints. For a 600 GeV Higgs mass, mixing angles between 0.001 and 
0.04 are excluded, where for larger tan/3 the upper limits are slightly less stringent. For a 



13 Cf [40] for a generic introduction, and 39 for recent work. 

14 We want to mention that larger heavy Higgs masses allow for Ai > for running up to arbitrary scales, cf. e.g. 
|36l |4"T] . However, the mass of the second Higgs boson is typically much above the LHC reach in the according 
setup. We thank O. Lebedev for useful discussions regarding this point. 
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Running couplings, constrain! from perturbativity ofX 3 Running couplings, constraint from vaccuwn stability 




5 10 15 20 10 20 30 40 50 

t t 

FIG. 3: Example for RGE running: LEFT: A3 becomes non- perturbative at t ~ 21 [/x run ~ O(10 7 GeV)], 
with sin a = 0.37, tan/3 = 0.43, RIGHT: Example for region where 4Ai A2 — A3 only varies marginally 
over large scale ranges, with sin a = 0.12, tan/3 = 0.04. 

Higgs mass of 1 TeV, this region is decreased to 0.02. Increasing sin a, the transition into 
the allowed region comes from an enhanced value of Ai at the low scale; in this case, the 
limiting value is again the perturbativity of Ai. Note that in parts of the parameter space 
A2,A3 only change marginally; in this case, there is a very fine interplay between the rise 
of the absolute values of A2, A3 and the rapid decrease of Ai, so including additional orders 
in the running might change these bounds, leading to a larger allowed region. An example 
for such a "slow-running" point is given in Figure [3} If we want to prevent this fine-tuning 
over large scales, we could e g allow for slightly negative values of 4Ai A2 — A§; opening 
up the condition such that 4A1A2 — A3 > —0.001 leads to (sina) m i n ~ 0.015(0.01) for 
a 600 GeV (1 TeV) Higgs mass. In priniciple, this area of parameter space would need a 
more detailed investigation. However, this region is phenomenologically difficult test, and 
the most important limits are indeed the ones from perturbativity on the maximal allowed 
mixing, so we will not investigate this in more detail in this work. 



Summary of RGE effects 

In this subsection, we will first summarize the results for a 600 GeV Higgs at a running scale 
corresponding to t = 37 and then discuss variations of the heavy Higgs mass and consequences 
when going to a higher scale. In Figure [4j we show the allowed parameter space for 600 GeV Higgs 
mass both at the low (t = 37) and the Planck scale. As discussed above, the largest constraints on 
large mixing angles are given by running of Ai and A2 for low/ high tan f3 regions respectively, while 
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FIG. 4: Limits for m# = 600 GeV at t — 37 (left) as well as Planck scale (right). We here consider 
| sin a | < 0.49; the experimental limit is given by |sina| < 0.23 

generally tan/3 > 0.37 (0.26) is excluded by A2 running at the low (high) scale. Additionally, small 
| sin a | values are generically excluded from requiring vaccuum stability. The minimal/ maximal 
values for |sincv| are ~ 0.035/0.3(0.1/0.2) at t = 37 (the Planck scale). Although the above 
discussion focuses on a Higgs mass of 600 GeV, the characteristics of the respective limits remain 
the same if the mass or the scale of the running are increased. We observe the following effects: 

• Raising the heavy Higgs mass while keeping the scale fixed leads to a reduction of the maximal 
allowed mixing angle, which stems from the perturbativity of Ai, as well as a decrease of 
the allowed maximal value of tan/3 from perturbativity of A2. However, on the other hand 
smaller mixings are still allowed. This is due to a larger A2 value at the EW scale, which 
prevents a fast decrease of 4Ai A2: this equally holds for larger A3. Even for negative j3\ x 
function values at low scales, the growth of A2, A3 can prevent Ai from becoming negative 15 . 
In general, the allowed region shrinks and equally moves to smaller mixing angles and tan f3 
values. The effects are displayed in Figure [5j where we compare the allowed parameter space 
at the Planck scale for a 600 GeV as well as 1 TeV heavy Higgs mass. 

• Raising the scale while keeping the Higgs mass fixed has similar effects: the maximal allowed 
mixing angle area is further restricted; generally, the allowed region is shrinking and moving 
to smaller minimally allowed tan j3 values, cf Figures [4] and [5} 



E.g., such a point is given by sina = 0.1, tan/3 = 0.05, which is excluded (allowed) for mjj = 600GeV(lTeV) 
by requiring stability up to the Planck scale. 
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allowed parameter space for different masses at Planck scale allowed parameter space for different scales 




sin a sin a 

FIG. 5: Limits at Planck scale for different to# (left) as well as limits for to# = 600 GeV at different 
scales (right). SM signal strength measurements have been neglected. 



allowed regions for varying Higgs masses at t=37 scale allowed regions for varying Higgs masses at t=37 scale 
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FIG. 6: Contour plots for different (fixed) m# at t=37. In the right plot, the third vacuum stability 
condition 4A1A2 > A§ has been neglected, leading to a slightly enhanced region for small |sina| values. 
SM signal strength measurement has been taken into account. 

In Figures [6] and [7j we present the results of our scans including all limits in terms of contour plots 
for the allowed areas at t = 37 as well as the Planck scale for = 600, 700, 800, 900, 1000 GeV, 
with numerical values summarized in Tables [I] and |TT{ As discussed above, the validity of the third 
vacuum stability condition, i.e. 

4 Ai A 2 - A3 > 

using NLO precision only might be questioned, so we equally present results where this is neglected. 
In general, this opens up the parameter space for even smaller mixing angles. We then take the 
constraints from vacuum stability following Ai running as a conservative lower limit. 
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allowed regions for varying Higgs masses at the Planck scale allowed regions for varying Higgs masses at the Planck scale 




sin a sin a 

FIG. 7: As Fig. [6j but for validity up to the Planck scale. 



m H 


sina| m i n ,37 


sin a maXj 37 


(tan^) maXj 37 


600 


0.035 


0.23 


0.37 


700 


0.030 


0.23 


0.31 


800 


0.024 


0.21 


0.28 


900 


0.019 


0.19 


0.25 


1000 


0.016 


0.17 


0.22 



TABLE I: Allowed parameter ranges for varying Higgs masses at the low (t = 37) scale; | sina| m i n max 
taken at tan/? = 0.15. x < 1 TeV fixes the lowest tan/3 value to 0.025. For run < 700 GeV, the maximal 
allowed mixing angle results from the measurement of the light Higgs signal strength. 

IV. TRANSLATION TO COLLIDER OBSERVABLES 

The parameter space presented in the last subsection translates into two different observables at 
colliders: 

• the generic suppression of the production of the heavy Higgs; this is given by sin 2 a, 

• suppression of the SM decay modes of the heavy Higgs. Here, we have to take into account 
that the additional mode 

H hh 

leading to a further reduction of the SM-like branching ratios. 
The total width of the heavy Higgs is then modified to 

Ttot = sin 2 aTsM, tot + ^H^hh- 
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777 zjt 


sin rv • ni i 
I "Imm, Planck 


sin rv\ m i 
I u | max, Planck 


(tan /?i ni i 
V tclAi M 'max, Planck 


Ron 


n i pi/i 




u.zo 


700 


0.086 


0.17 


0.22 


800 


0.074 


0.15 


0.20 


900 


0.064 


0.13 


0.17 


1000 


0.055 


0.12 


0.15 



TABLE II: Allowed parameter ranges for varying Higgs masses at the low Planck scale; | sina| m i nimax taken 
at tan/3 = 0.08 and at sin a < 0. As before, the minimal value of tan/3 is determined by the scan range. 

Following the observables tested by the experiments, we therefore consider 

with the global SM-scaling factor defined as 

o"bsm x BRbsm sin 4 ar to t,SM 

k = = — 

cjsm x BRsm r to t 
In analogy with the above definition, we also introduce a scaling factor k' which parametrizes the 
H — )• hh decay: 

/ obsm x BRbsm sin 2 aT H ^ hh 
« = = r ' 

CSM I tot 

where k + k! = sin 2 a. 

For a better understanding of the effect of the constraints on the (r, k) parameter space, we 
first investigate the H —> hh branching ratio. Figure [8] shows the constraints in this decay 
width from RGE running to the scale defined by t = 37. We see that different regions are 
excluded, depending on the sign of sin a, where the biggest effects stem from perturbativity of 
A2 • As discussed in Section [XTJ, the H — >• hh squared coupling is approximately even under 
a sign change of sin a for small tan/3 values; if tan/3 increases, it is larger for positive sin a 
values. In addition, for positive (negative) sin a values, the absolute value of the coupling 
decreases (increases) for decreasing tan/3. Taking this into account, the exclusion bounds from 
perturbativity of the couplings green plotted in Figure [8] are clear to interpret: for positive 
values of sin a, larger tan/3 values lead to larger decay widths, and therefore large values of 
^H^hh are here excluded from A2 perturbativity, while smaller ones are affected by Ai limits, in 
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Partial width bounds 









H A, 2 > 4 ji 








allowed 




J 



T 15 




0.4 -0.3 -0.2 -0.1 0.1 0.2 0.3 0.4 



FIG. 8: Exclusion of partial widths from perturbativity at the low scale for | sina| < 0.49. 

accordance with the limits in Figure [3j For negative sin a values, the roles of Ai and A2 are inter- 
changed, while smaller values of Th -+hh are achieved due to smaller absolute values of the coupling. 



Figures [9] and 10 then show how the cuts translate on the allowed parameter space in the (r, k) 
plane. The most important result is that the limitation of the angle comes with a maximally 
allowed total width Ttot < 14 GeV (reducing to < 6 GeV if we require perturbativity up to the 
Planck scale) for a Higgs mass of 600 GeV ; this should be compared with the SM-like Higgs Boson 
width of ~ 100 GeV. For a 1 TeV Higgs, the maximal values are 25 GeV for perturbativity at the 
low scale and 12 GeV for perturbativity at the Planck scale, respectively. Vacuum stability cuts out 
lower regions of the mixing angles/ Tu-^hh't here, it is important to note that the minimally allowed 
widths are in the sub-GeV range. The limits from this on the (r, k) plane are more pronounced 
for higher scales and lead to minimal values of r m i n = 0.79 GeV (1.4 GeV), ^min = 0.004(0.001) 
for validity up to the Planck scale for a 600 GeV (1 TeV) heavy Higgs. 



Finally, we remark that, if the limits from RGE are not considered, all other constraints are 
much less stringent; only the upper limit on tan /3 from perturbative unitarity cuts out a small 
region in the large positive sin a/ large T and large negative sin a/ low V region, similarly to the 
constraints that we obtain from perturbativity requirements of A2. 
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Constraints from perturbativity of couplings Constraints from perturbativity of couplings 




r r 

FIG. 9: Constraints on the (T, k) plane from perturbative running of the couplings (to# = 600 GeV) to 
t = 37. LEFT: Constraints for | sina| < 0.49; RIGHT: Constraints when the light Higgs signal strength is 
taken into account. 



allowed scale factor and total width, t=37 allowed scale factor and total width, Planck scale 




FIG. 10: Exclusion bounds in (r, k) plane from perturbativity, different heavy Higgs masses, for RGE 
running up to t = 37 (left) as well as the Planck scale (right). 

We find that the maximally allowed values for k are roughly (0.04; 0.04; 0.04; 0.03; 0.025) for 
n^H = (600; 700; 800; 900; 1000) GeV. Concerning collider searches, the best prospect is therefore 
the search of a relatively light Higgs Boson at mn = 600 GeV, which would lead to 0.013 (0.093) pb 
at a 7(14), TeV LHC in the gluon fusion and 2 x 10 -3 (0.016) pb for the vector Boson fusion 
channel 16 . However, also note that the widths for all masses are < 25 GeV, which might allow 
for new search strategies at such masses for narrow scalar resonances 17 . The maximal value of «/, 
on the other hand, is 0.013 for a 600 GeV Higgs mass, which would lead to a total cross section 

16 Production cross sections have been taken from |42j . 

17 Present studies usually assume quite broad Higgses in this mass range, following the SM Higgs searches. This 
assumption is not consistent with our scenario. 
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Constraints from perturbativity of couplings 



allowed scale factor and total width, t=37 





FIG. 11: Exclusion bounds in (r, k') plane from RGEs, t = 37. LEFT: effect of perturbativity of different 
couplings. As before, the strongest constraints stem from the requirement of perturbativily of A2. RIGHT: 
Results when varying m#. The maximal value of n' is in the % range for run = 600 GeV. With increasing 



of 4 x 10 3 pb(0.03pb) from gluon gluon fusion production at a 7 TeV(14 TeV) LHC for the 
additional channel H — )• hh. Other allowed values of (k' ,T) for different Higgs masses at the low 



(t=37) scale can be obtained from Figure 11 If, as briefly mentioned in Section HID, we relax 
the requirements of both perturbative unitarity as well as electroweak symmetry breaking at high 
scales, this basically opens up the parameter space for smaller positive mixing angles, effectively 
leading to lower minimal values of n, T. The effects are negligible for the low scale; for the Planck 
scale, the minimal allowed width is decreased to ~ 0.5 GeV for nearly all masses considered here. 
However, as the small mixing range will be hard to detect at colliders, there is no visible impact 
from this on the above discussion of collider observables. 



Variation of input parameters 

Finally, we comment on how a variation of SM-input parameters affects our results. The running 
of the SM Higgs coupling is known to be sensitive, especially 18 to the strong coupling a s as well 
as the top Yukawa coupling yt, so we investigate the robustness of our results under variations of 
these parameters at the low scale. We here consider 

a s (m z ) = 0.1184 ± 0.0007, y t (m t ) = 0.93587 ± 0.002 GeV, m t = 173.1 ± 0.7 GeV, 



Cfe.g. [H]. 
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where, following [39], we combine the above error in the top sector to 19 

y t (m t ) = 0.93587 ± 0.006 GeV 

For the results presented so far, we have used the central values above. We comment directly on 
the effects of values for these are the observables of main experimental interest. 

• lowering (raising) a s , while keeping yt fixed: 

in this case, the allowed minimal mixing angle allowed from vacuum stability is marginally 
increased (decreased). A changed value of a s mainly influences the running of yt, which in 
turn leads to a faster (slower) descent of Ai, mainly affecting limits from the third vaccum 
stability condition. Changes are however in the % regime. The upper limits of sin a, which 
constitute the main restrictions in the (r, k) plane, are not affected. 

• keeping a s fixed, while raising (lowering) yt: 

for a higher (lower) yt, a larger (smaller) region of small mixing angles is excluded, again due 
to the faster (slower) descent of Ax- Maximally allowed values as well as large tan/3 limits 
are not changed. However, using yt(mt) = 0.92987 as input value shifts the breakdown of 
the SM-case (sin a = 0) running of the coupling by approximately an order of magnitude 
to t ~ 42.6. Limits from perturbativity of the couplings still persist. 

V. CONCLUSIONS 

In this work, we have investigated the theoretical and experimental limits of the parameter 
space of a pure singlet extension of the SM Higgs sector without contributions from a hidden 
sector, where the heavy Higgs lies in a mass range of 600 GeV to 1 TeV. We found that, after 
the light Higgs coupling strength measurements from the LHC experiments have been taken into 
account, additional strong limits stem from perturbativity of the couplings as well as vacuum 
stability, following from the /3-functions of the theory. Even for a relatively low breakdown scale 
O(10 10 GeV), the running of the heavy and light Higgs self-couplings severely restricts the allowed 
parameter space. We have translated this into observables which are currently tested by the LHC 
experiments, i.e. a global rescaling factor k for SM-like decay modes of the model, as well as total 
width r of the new scalar. In the heavy Higgs mass range considered in this work, the light Higgs 



Note that we only want to estimate the effects of deviations from the central values; a more accurate error 
determination should be embedded in a higher order investigation of the Higgs singlet extension. 
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Boson signal strength restricts k to 0.04 for ran ^$ 700 GeV (at the low scale), while for higher 
masses additional constraints arise from the running of the couplings . Hence, the searches for 
such a Boson at the 7/8TeV with a relatively low luminosity are surely challenging. However, 
on the upside we found that the total width of the new scalars is usually quite suppressed with 
respect to SM Higgses of such masses, with widths lying in the 1 — 25 GeV range (they are always 
< 0.02 ran)- In addition, we have introduced a second scaling parameter k' which parametrizes 
the additional decay H — > hh. We found that maximal values of this parameter are in the % range. 

In our work, we have neglected additional contributions in the /3-functions which might modify 
the runnings and eventually enhance the parameter space in the tan /3 > 0.2 region stemming from 
the running of A2- These contributions, which could originate from the hidden sector, would have 
to be large and negative, cancelling the rapid rising of the couplings which leads to the exclusion of 
experimentally interesting regions with scaling factors k being limited by the mixing angle alone. 
Scenarios with larger n' values which parametrizes Tn^hh are equally suppressed by the running 
of A2- We plan to investigate such options and the corresponding phenomenological implications 
in future work. 
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Appendix A: Discussion of Higgs potential and vacuum stability conditions 



In this section, we briefly guide the reader through the steps from Eq ^ to Eq ([3]), using the 
definition of the scalar fields given in Eq ([5]). Since the minimisation procedure is not affected by 
the choice of the gauge, it is not restrictive to define the two VEVs in the following way: 

v_\, {x) = ^> (A1) 

VV27 

with v and x real and non-negative. 

Then, the search for extrema of V is made by means of the following differential set of equations: 



dV 
dv 
dV 
k dx 



(v,x) 
(v,x) 



X 



2 , \ 2 , -^3 2 

-m + Aiv H x 



2 , \ 2 , A3 2 

-[l + \2X + —V 



(A2) 



The physically interesting solutions are the ones obtained for v, x > 0: 

v 2 



X 



A2m 2 




A1A2 


A 3 
4 


Ai^ 2 - 


- ^m 2 


A1A2 


\ 2 

4 



(A3) 



(A4) 



Since the denominator in equations (A3)-(A4) is always positive (assuming that the potential 



is well-defined), it follows that the numerators are forced to be positive in order to guarantee a 
positive-definite non- vanishing solution for v and x. 

In order to identify the extrema, we need to evaluate the Hessian matrix: 



/ d 2 v d 2 v\ 



H(v, x) 



dv 2 

d 2 v 



dvdx 

d 2 v 



' 2\ lV 2 



\-zvx 



\ 



\%vx 2X2X' 



(A5) 



V dvdx dx 2 I 

From this equation, it is straightforward to verify that the solutions are minima if and only if 
equations ^ are satisfied. 

To compute the scalar masses, one must expand the potential in equation ^ around the minima 
found in equations (A3)-(A4). Then, Eqns. follow immediately. 



Appendix B: Perturbative unitarity 

In this section we want to briefly explain the techniques that we used in order to obtain bounds 
from perturbative unitarity, firstly described in detail by [25J. Evaluating the tree-level scattering 
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amplitude of longitudinally polarised vector bosons one finds that the latter grows with the energy 
of the process, eventually violating unitarity, unless one includes some other (model dependent) 
interactions. According to the equivalence theorem, the amplitude of any process with external 
longitudinal vector bosons Vl (V = W^, Z) can be substituted each one of them with the related 
Goldstone bosons v = w^, z [27] for energies much larger than the vector Boson mass. 

Given a tree-level scattering amplitude between two spin-0 particles, M(s,9), where 9 is the 
scattering (polar) angle, we know that the partial wave amplitude with angular momentum J is 
given by 

i r 1 

aj = ^r d(cos 9) Pj (cos 6) M(s, 6), (Bl) 
62n J_ 1 

where Pj are Legendre polynomials. It has been proven (see [26]) that, in order to preserve 
unitarity, each partial wave must be bounded by the condition 

|Re(aj( S ))| < i (B2) 

As discussed previously, in the high energy limit, y/s — > oo, only the ao partial wave amplitude 
does not vanish; therefore, we here present all ao's in the high energy and small gauge coupling 
limit (i.e. y/s — > oo and e — > 0, respectively): 



^ = 3(rn 2 h + rn 2 H + (rn 2 - m 2 H )co S (2 a )) 

64ttv z 



m\ + m 2 H + (m^ — m? H ) cos (2 a) 



cos a , . 2 2 
a (zz — >■ hh) = — — — \{om h + m H )xcosa+ 

1207TV Z X 

+ (m\ — m 2 H ) (x cos(3 a) — 4v sin 3 a)) (B5) 

cosasina((m| + m 2 H )x + (m\ — m 2 H )(xcos (2a) + -usin(2a))) 
a (zz -+ hH) = ^- 7 =- 5 - (B6) 

sin a , . 9 2 s 
ao(zz^-HH) = — — ^— (2(— m h + m H )v cos a+ 
G4irv z x 

+ x(m 2 l + m 2 H + (rr? h — m 2 H ) cos (2 a)) sin a) (B7) 
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, . m( + mir + (mi — mir) cos (2 a) , n . 

a (w + w~ -> w + w~) = — * 2- v h H ' v - B8 

lo7rtr 

ao(w + w~ —> hh) = -= {(3ml + m 2 H )x cos a+ 

64\/2irv 2 x 

+ (m 2 h - m 2 H ) (xcos (3 a) - 4v sin 3 a)) (B9) 

( + - >im cos a sin a((m 2 + m^)x + (m| - m^)(a cos (2a) + v sin(2a))) 
a (w u; -+hH) = 32~^x ( ' 

a (w+w- ^ HH) = —^^—(2(-ml + m 2 H )vcos 3 a+ 

+ x(ml + m 2 H + (ml — m 2 H ) cos (2 a)) sin a) (Bll) 



a (hh^hh) = Yo2 ±^(6(5m 2 h + m 2 H )(v 2 + x 2 ) + 

- 3(15ml + m 2 H )(v — x)(v + x) cos (2a) + 
+ 6(3m 2 - m#) (t> 2 + x 2 ) cos(4a) + 

- 3(m 2 h -m 2 H )((v-x)(v + x)cos(6a) + 8vxsm 3 (2a))) (B12) 

3 cos a sin a . 2 2 

ao(nh^/-hH) = — — — — -(xcosa + v sma)((3m h + m^jxcosa + 



64t/2itv' v 



+ (ml — m 2 H )x cos (3 a) — (3ml + m 2 H )v sin a + 

+ (m^ — m 2 H )v sin(3a)) (B13) 
a (hh ^ HH) = (Q(-m 2 h + m 2 H )vx cos(4 a) + 

+ 6(ml + m 2 H ) (v 2 + x 2 ) sin(2a) + 

- (m 2 h - m 2 H )(2vx + 3(v - x)(v + x)sin(4a))) (B14) 



a (hH ^ hH) = ^^^^ (6(-ml + m 2 H )vx cos(Aa) + 
+ 6(m 2 h + m 2 H ) (v 2 + x 2 ) sin(2 a) + 

— (ml — m 2 H )(2v x + 3(v — x)(v + x) sin(4a))) (B15) 

/, tt tttt\ 3 cos a sin a . 2 2 \ 

ao(nH ^ HH) = = — — — (vcosa — a;sina)((m h + 3m i? )vcosa + 



64\/27rr- r 



+ (— m fe + m H )vcos (3 a) + 
+ 2x(ml + m#- + (m^ — m 2 ^) cos (2 a)) sin a) (B16) 
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aa(HH — »■ HH) = —±^(6( m 2 +5m 2 H )(v 2 + x 2 ) + 
+ 3(m^ + 15m%)(v — x){v + x) cos (2 a) + 

— 6(m| — 3m^) (v 2 + x 2 ) cos(4 a) + 

- 3(ml-m 2 H )((v-x)(v + x)cos(6a) + 8vxsm 3 (2a))) (B17) 

Appendix C: Analytic solution for SM gauge coupling RGEs 

In the SM, all one-loop RGEs for gauge couplings are of the form 

dx o 
— = a x . 
dt 

The exact analytic solution for this equation is given by 

*(*)=i v (CI) 

v ; l-ax(t = t ) (*-<o) 



where for t = log ( j^- ) we have 



rc f 



t - t = 2 log 



For positive values of a, the coupling reaches the Landau pole when the denominator in Eq. (CI) 
goes to 0; for negative values, x — > for t — > oo. 

Now we turn to the Yukawa coupling terms. This generic equation has the form 

dx , o 

— = ax + ox 
dt 

with the solution 

x , t) = Vamto)e a{t ~ to] 

y/l-bet'tt-^C'ito)' 

2 

with C'(to) = t h 2 where x(t = to) = defines the initial value. In case of the top Yukawa 



coupling, we have 



167r 2 a = -Ag 2 s --g 2 --a' 2 ] 16vr 2 6= - A . 

8 24 4 



However, taking the time dependence of the SM gauge couplings into account, the above solution 

ha 



needs to be modified such that a (t — to) is replaced by f. a(t') dt'. Although this is still feasible 



at one loop, we chose to solve the RGE of the top Yukawa coupling numerically 
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Cf also [44] for an all-analytic solution to the first order RGEs. 
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